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Abstract 

The extended Hamilton's Priciple and other methods proposed to handle non-holonomic 
constraints are considered. They dont agree with each other. By looking at its consis- 
tency with D'Alembert principle for linear non-holonomic constraints, it was claimed 
in earlier papers that the direct extension of hamilton's principle is incorrect. Non- 
holonomic Constraints, linear in velocities were considered for this purpose. This 
paper analyzes these claims, and shows that they are incorrect. And hence it shows 
that it is theoretically impossible to judge any attempt on non-holonomic constraints 
to be wrong, as long as they are consistent with the D'Alembertian for holonomic 
constraints. 



I. Introduction 

The problem of mechanics with non-holonomic constraints was considered re- 
cently in [21 [2] and [5]. pQ and [3], appear to be among the earliest to propose a 
direct extension of the hamilton's principle to the non-holonomic regime. The 
same extension was adopted in [15]. However, later papers ([3] [HIE] to name a 
few) attempted to show that such an extension is incorrect. [3J appears to be 
very conclusive regarding this issue. It has ruled out the extension of hamil- 
ton's principle usimg a theoritical argument. The present paper shows that such 
theoritical arguments cannot rule out the extension. 

This paper discusses the extension of hamilton's principle along with other 
attempts to tackle non-holonomic constraints, and carefully examines the argu- 
ment used to rule out the direct use of hamilton's principle for non-holonomic 
constraints. Finally, it points out the reason why the argument is incorrect. 
That slightly alters the conclusions made in [3J regarding the scope of the 
D'Alembert's principle. These conclusions show that the two ways of tackling 
non-holonomic problems are uncomparable. 

The paper first briefly discusses the classical theory of mechanics, which can 
be found in any classical text on mechanics [12J. However, the following section- 

II outlines the theory in a structure convenient for the present purpose. Section 

III analyzes the various proposals to tackle non-holonomic constraints. Section 

IV analyzes the argument given against a direct generalization of the Hamilton's 
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principle, and points out the fallacy in it. Section V summarizes its implications. 
The mathematical tools required, namely, the variational calculus is put at the 
end, in the appendix. Einstein's notation is used for summation throughout. 
Unless explicitly mentioned, repeated indices mean summation throughout. 

II. Holonomic Constraints 

1. Newton's Laws 

For a system of particles Newton's law can be written under a Cartesian system 
as: (indices not summed) 

Where Xi are the total forces. Newton's law says, acceleration is caused by force 
and only force. Hence, within the framework of newtons laws, nothing other 
than a force can influence the dynamics of a system. Any constraint which 
restricts the motion of a system should hence be replaced by an equivalent 
'constraint force' which brings in the same effect as the constraint, so that it 
can be studied using Newton's laws. The equation becomes mii\ = Fi + Qj, 
where Fi are the applied forces andQi are the constraint forces. An additional 
postulate is required to construct the appropriate constraint forces. 

2. Constraint Force Postulate (D'Alembert's Principle): 

For a constraint of the form ^{x,t) — 0, which means, that the system is con- 
fined to a time varying hyper surface in the configuration space, the constraint 
forces Qi which have the same effect as that of the constraint can be written as, 

Q. = A(t)Aj? (i) 

where X(t) is an undetermined function of time. This formulation of the pos- 
tulate will be shown to be equivalent to the standard D'Alembert's principle 
stated in terms of virtual work, defined in most standard texts on mechanics ( 
|12j . for example). 

A 'virtual displacement' is a local variation of x(t) consistent with the con- 
straint, i.e., a local variation {5i(t)} such that, J? Xi $i = , assuming that not 
all jF Xi are zero, at a time t. This is a simpler form of the definition given by 
[16] , It says "virtual displacement is a vector tangential to the constraint man- 
ifold" The definition of a virtual displacement is found all texts on mechanics. 
However, most of them are unclear and ambiguous. Hence we refer to [15]. As 
it turns out, the major result pointed out by this paper banks on this defini- 
tion, hence this extra care on the definition. D'Alembert's principle says the 
virtual work done by the constraint forces vanishes|12|. i.e., QiSi — 0, at a time 
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t for all {Si(t}} such that ,^ Xi 5i — Or, Qi = X.^ Xi where the constant A 
might not be the same at a different time t'. This can be in general written as 
Qi = \(t),^ Xi . Hence equation Q] is equivalent to D'Alembert's principle. The 
equations of motion can now be written as 

rrnxi = F t + \{t)& Xi (2) 

where F$ are the applied forces. 



3. Variational Formulation 

For a system where the applied forces are conservative, the equation of motion 
[5] supports a variational formulation. If jSf is the Lagrangian, the equations [5] 
can be written as (|12|): 

j t ^ Xz -^ Xl =\{t)& Xt ) (3) 

Since & is independent of the velocities, \{t)^ Xi = -{^|gf:(A(t) t F)-gf-(A(t) I ^')} 
. Hence the equations can be written as -^-£-(3? + \(t)^) - ^f-(-S? + A(i)^") = 
0. Which can be compared with equation [21 in the appendix. These equations 
suggest the Hamilton's principle of stationary action for this category of con- 
straints. The equation l21l along with initial and final condition (i.e, the boundary 
value problem) is an extremization of the action. However, the corresponding 
initial value problem (i.e, the same differential equation with initial position 
and velocities known, which is the case in all physical situations) is not strictly 
equivalent to extremization of the action. However, the Hamilton's principle 
can be formulated with the obtained boundary values as 'action is stationary 
for fixed boundaries'. 



III. Non Holonomic Constraints 

The D'Alemberts virtual work postulate stated in the previous section fails 
to work with nonholonomic constraints. A new constraint force postulate is 
required to handle non holonomic constraints. Various postulates are available 
in literature [2]. The major candidates are: 



1.) Gauss-Gibbs Principle 

The Gauss-Gibbs principle as described in [Ml [6], says that the first order 
variation of the quantity, C, vanishes under allowed Gaussian variations. 

C = lmi(xi - — ) 2 (4) 
2 n mi v ' 
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Variations {<$,(£)} for which Si and Si vanish at the concerned value of t are 
called Gaussian variations. Under such variations, the first order variation of C 
is given by: 

5C = (m t x'i - Xi)5i (5) 

This quantity vanishes for all Gaussian variations consistent with constraints. 
The consistency condition translates as: 

&(x + 5,x + 5,t)- &(x, x, t) = & Xi 5i + & Xi S\ = (6) 

Which is trivially true at the concerned time t. And hence we turn to its first 
order variation in time. 

6&(t + St) - S&(t) = {./, S, + (J^ + & Xi )S\ + & Xi Si}8t = (7) 

That puts the condition & ' Xi Si = on the Si (Second order time variation needs 
to be invoked for holonomic constraints). The principle is now equivalent to 
newton's laws with a force postulate: 

Qi = \{t)^xi (nonholonomic) (8) 

Qi = \(t)^ Xi (holonomic) (9) 

A similar result is obtained in [3]. 

2.) Jordain Principle 

The fundamental equation of the Jourdain principle is([71 114j) 

{m lX \ - Xi)Si = (10) 

Where Si(t) are a 'Jourdain' variations which vanish at t, while their deriva- 
tives dont. Such variations are called 'virtual velocities'. The constraint force 
postulate used is same as those in the Gauss's principle (JSEJ) 



3.) Extended Hamilton's principle: Vakonomic Mechanics 

Extending Hamilton's principle of stationary action was proposed in [TJ[H], later 
criticized heavily ( [2J [31 E] for example.) Nevertheless, it is equivalent to assum- 
ing a constraint force given by 

Qi = A(i){J^-^J + J^^ (11) 

This was termed as 'vakonomic mechanics' (Variational Axiomatic Kind) by 
|17| since the basic axiom here is the Hamilton's principle, which is variational 
in nature. 
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We now have two possible postulations of the the constraint force: IH0and[TTl 
We have one well accepted formulationQjD'Alembert's principle) for holonomic 
constraints. Both 00 an d [TTJ agree with [TJ for the special case of holonomic 
constraints. But, they don't at proper nonholonomic constraints. One curious 
and simple set of constraints which might be of help in deciding the correctness 
of00and[TT]is the linear nonholonomic constraints. This has been used in 
to rule out [TTJ However, there seems to be a subtle issue with that and hence 
we wish to reconsider it. 

IV. Linear Non-Holonomic Constraints 

Linear nonholonomic constraints are non integrable equations of the form, 

a l x i + a t = (12) 

Where at and at are functions independent of ii. This category of constraints 
has been recently considered in [3] to discard [TTJ [3] has concluded that all 
nonholonomic constraints are beyond the scope of Hamilton's principle; the 
general nonholonomic constraints are beyond the scope of D'Alemberts princi- 
ple, while the linear non holonomic constraints do remain within the scope of 
D'Alembert's principle. And they agree with what 00 says; and [TTJ is inconsis- 
tent with it, hence it is incorrect. This has been done in other literature as well. 
[3J examines the reason behind it's conclusions. And it says, the reason why [TTJ 
fails at a linear non-holonomic is that the allowed variations are not consistent 
with the constraints. By allowed variations, we mean the variations constrained 
by 

aA = (13) 

which is indeed obtained from the constraint equation, but not directly stating 
that the variations are consistent with the constraint. But, the constraint force 
given by [I3J is Qi — A(t)aj. Which is same as what is obtained from 00 as 
mentioned above. However, 00 obtaines the constraint force through a different 
route, not from 1131 Since equation [13J is mysterious, we take a look at how to 
arrive at [T3J from [T2T 

It is widely mentioned that [T3] is obtained from [T^] using D'Alembert's prin- 
ciple of virtual work([U [31 HOI US HH US)- They writeDJas 

aidxi + atdt — (14) 

And then replace dxi by 8i and dt by 0. But the trouble is, this cannot be done, 
if we define a D'Alembertian displacement, Si(t) as a variation in Xi(t). dxi is 
a differential motion over a time dt, and is given by dxi — iidt. It vanishes if 
dt is set to zero. This definition of virtual displacement is used in |16| . A true 
expression of the constraints on the variation is, 



aiSi + a iXj XiSj + a tX] S 3 = 



(15) 
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D'Alembertian displacements don't say anything about the derivatives of the 
displacements. They are variations in position, by definition. As shown earlier, 
the same constraint force can be obtained by a different route using Gaus- 
sian principle or the Jourdain principle. But the force there is expressed as 
Qi = X(t)^xi for non holonomic constraints, where & = ciiXi + a t . More- 
over, Jourdain displacements are virtual velocities, not virtual displacements. 
Hence it is impossible to deduce the constraint force Qi — \{£)&xi directly 
from D'Alembert's virtual displacement principle, even for the case of linear 
non holonomic constraints. 

Hence, equation [T3J which is the basis of all support to 00 an d all criticism 
on[TT]is incorrect. It means, the linear non holonomic too, is outside the scope 
of D'Alembert's principle, and hence is no way of deciding between IHUl and [TO 
There is no theoretical way to decide between 00 an d HU Neither of them 
are intrinsically incorrect. However, 00 seems to be accepted in most practical 
examples. 

V. Concluding Remarks 

To summarize, there are two approaches to mechanics. One based on a con- 
straint force postulate, namely the Jourdain/Gaussian principle(that covers the 
D'Alembertian also) and the other based on a variational principle, namely, the 
Hamilton's principle. Mechanics based on Hamilton's principle is called Vako- 
nomic mechanics. More theoretical aspects of the Vakonomic mechanics can be 
found in|17]. The D'Alembert's principle is the most intuitive among all. It uses 
the intuitive concept of a constraint force acting instantaneously normal to the 
constraint surface, and hence forcing the particle to the surface. And hence it 
is most reliable. Earlier papers have attempted to compare the two approaches 
by checking their consistency with the D'Alembert's principle. In this paper, 
we showed that such a comparison is impossible. However, with an application 
point of view, the constraint approach is found to be more oftenly used. 

VI. Appendix: Calculus Of Variation 

A detailed account of calculus of variation can be found any standard text on 
the subiectfTgl [T51 [2U] for example). This section contains a brief outline of 
calculus of variation in a convenient format, including all the results relevant 
for this paper. This whole section refers to |18[ 119) . 

Calculus is the study of local variations. For a functional Jz?(f,f',t) of a 
function f(t), and its integral defined as 




(16) 
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where /' is the time derivative of /. The local variation in /(/) for a local 
variation of fi(t) from fi(t) to fi(t)+5i(t) is given by 

ti t 2 
I(f+S)-I(f) = J (J? fi 6 i+ J? f ,6l)dt+ J (^f^+^f^S'j+^f^S^dt 
ti ti 

(17) 

up to the second order. Under the additional condition that /i(ii)and fifo) 
are fixed during the variation i.e., Si(ti) = Sifa) = 0, the first order variation 
in 1(f) can be written as, 

t 2 

Hf + Hf) = J {# fl j t (^f>Wt (18) 
ti 



1.) Unconstrained Extremization 

Extremization of 1(f) would mean that finding a f(t) where its first order vari- 
ation vanishes, i.e, J t ^{^'f i — -^(Jz? f>)}Sidt = Ofor all functions Si. That is 
equivalent to the n differential equations -5?/* — (^ /; ) = with the boundary 
conditions given by /j(£i)and /i(^) 

2(a) Constrained Extremization. 

We first consider a constraint of the following type: J t * 2 J?(/, /', t)dt = where, 
&(f,f',t) is a functional. This problem is now equivalent finding functions 
{fi(t)} such that, 

(i) They satisfy the constraint J?(f,f',t) 

(ii) First order variation of 1(f) vanishes under local variations in{/i(i)} con- 
sistent with the constraint &(f, f',t) and the boundary conditions on f(t) 

(iii) { /i (t ) }satisfy the boundary conditions. 

Let {8i(t)}be a set of n functions and {e^} be a set of n real numbers. 
Condition (ii) can be written as: 

t 2 

J{Jzf fi -±(& f ,)}S i e i dt = (19) 
ti 

for all {ej G S where S = {{ej : e, f . - ±(3? f ,)}Sidt = 0}. For all sets 

{5i(t}}. Essentially, S is the set of all variations consistent with the constraints. 
Such a definition of S will not work if «^/ 4 — ^ (& f) — whenever & = for all 
i. The second order variation should be invoked in such cases. For a particular 
set {Si}, let on = f£{# fi - f t {^ r )}Sidt and ft = f£{& fi - f t (& f>)}Sidt. S 
is now the set of n- vectors {e^} orthogonal to the vector (3. Hence the condition 
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says, a,i — A/3j for a real number A. A simple argument can establish that the 
number A is independent of the set {Si}. Hence, we have the final equations: 
ai = X/3i for all sets {Si}. Or, equivalently, 



A is the Lagrange multiplier for this case. This is easily generalized to several 
constraints. 

2(b) Constrained Extremization: constraints of second type. 

Now we consider constraints of the form &(f, f',t) = 0. These constraints can 
be equivalently written as 

f£k(t)&(f,f',t)dt = for all funct ions k(t), or equivalently, for k(t) = 
kj(t), where {kj(t)} form a basis for functions on [ii,^] Now, the constraint 
has reduced to the previous category of constraints. The equations may now be 
written as: 



where X(t) is the Lagrange multiplier for this case. 
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